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Summary 

This report presents the conceptual development of a life extending control system where the objective is to achieve high 
performance and structural durability of the plant. A life extending controller is designed for a reusable rocket engine via damage 
mitigation in both the fuel and oxidizer turbines while achieving high performance for transient responses of the combustion chamber 
pressure and the raixtura ratio. This design approach makes use of a combination of linear and nonlinear controller synthesis 
techniques and also allows adaptation of the life extending controller module to augment a conventional performance controller of 
a rocket engine. The nonlinear aspect of the design is achieved using nonlinear parameter optimization of a prescribed control 
structure. 


1. Introduction 

Systems with high performance requirements and high power densities such as the Space Shuttle Main Engine (SSME), 
hypersonic propulsion engines, and gas turbine engines often have a small number of critical components that operate close to 
mechanical design limits. The critical components are also indicators of the effective lifetime of the entire system. These components 
often experience maximum stress conditions during transients, and it is during such transients that large decrements in the component 
life are experienced. Possible damage modes include spalling, creep, corrosion, and fatigue. Simply minimizing stress levels is not 
always a solution to these problems because it will typically result in an excessive loss of dynamic performance. 

The effect of thermal transient loading on the Space Shuttle Main Engine (SSME) turbine blades during startup and shutdown 
is a typical example of the above scenario. It was against the backdrop of the durability problems of the SSME that the concept of 
life extending control (LEC) has evolved. The fundamental concept of life extending control is to control the rates of change and 
operating domains of some performance variables to minimize damage (or damage rates) of the critical components while 
simultaneously maximizing the dynamic performance of the plant. While the life extending control technology was developed 
initially for the SSME, it has broad application to many situations where both rapid response through transients and long life are 
required. 

The fundamental concept of life extending control has been forwarded by Lorenzo and Merrill (1991a and 1991b). The 
following basic approaches have been considered: (1 ) implicit life extending controls, which use current cycle based damage laws, 
and (2) direct life extending controls which assume the development of a continuous form of damage law. The availability of a 
continuum damage model allows a more straightforward development of the life extending control concept and, hence, a simpler 
implementation. In the life extending control implementations that have been considered, in addition to the plant and the performance 
controller, a structural estimator must be addended which provides the stress, strain, and temperature states of the critical 
components. These are used by an appropriate continuum damage model which in turn provides estimates of the current damage 
rates for the damage controller. A continuum fatigue damage model based on the local stress method has been developed by Lorenzo 
(1994). 

Ray et al. (1994a and 1994b) have shown that, in an open-loop setting, it is possible to reduce the fatigue damage rate and 
accumulation in the turbine blades of a reusable rocket engine (e.g., the space shuttle main engine) with little sacrifice in plant 
performance. Their damage reduction procedure, however, is based on an extensive off-line optimization and does not take 
advantage of on-line damage predictions or measurements. Also, the resulting feedforward signal is optimized for a particular set 
of initial conditions and a maneuver which must be specified a priori. This method may not be applicable to maneuvers and/or initial 
conditions not used in the optimization procedure. Dai and Ray ( 1 996) applied the same procedure to creep damage in the main thrust 
chamber wall of the same rocket engine. 
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Kallappa ct al. (1997) find an optimal open-loop control sequence which mitigates the creep and fatigue damage being 
i n the main steam header of a fossil-fueled powerplant. The feedforward control is augmented with a feedback 
controller to provide robustness to the system. 

Tangirala (19%) applies a combined feedforward/feedback controller synthesis methodology to a laboratory testbed. The 
synthesis procedure involves finding an optimal open-loop input sequence and augmenting it with a damage-mitigating output 
feedback controller. In a similar manner. Holmes, Tangirala, and Ray (1997) present a procedure for designing output feedback 
damage-mitigating controllers for a reusable rocket engine. 

Holmes and Ray (1997) use a fuzzy controller to mitigate damage in the turbine blades of a reusable rocket engine. A similar 
procedure is applied to a fossil -fueled powerplant in Holmes (1997). 

The approach taken in this report is to separate the design of the performance controller and the damage controller. That is, an 
aggressive performance controller is designed using standard (linear) techniques to achieve a high level of dynamic response for 
the plant, here a reusable rocket engine. Following this, the structure for the life extending control is added as an outer loop. It consists 
of a structural estimator, followed by a continuum damage model, in turn followed by a linear controller. The parameters of the linear 
control structure are then determined using nonlinear parameter optimization. The expected benefits of such an approach are the 
following: ( 1 ) the process is straightforward to apply, (2) it eliminates the necessity for determining the optimal open-loop response, 
and (3) it should be applicable to any form of input command desired. 

This report is organized into eight sections including the Introduction. In Section 2 a high level view of the life extending control 
system is presented. A description of the reusable rocket engine used in this study is given in Section 3. The damage model used 
to calculate fatigue damage in the turbine blades is presented in Section 4. Section 5 contains details on the design of the performance 
controller, and Section 6 describes the damage controller design procedure. Computer simulation results are presented in Section 7 , 
and the summary of the research and conclusions are given in Section 8. 


2. The Life Extending Control System 

The Life Extending Control (LEC) system functions as an integral part of the primary control loop, sometimes within the 
hierarchical structure of the intelligent control system (Lorenzo and Merrill, 199 la and 1991b). The focus of this section is on the 
fundamental issue of formulating a control structure for the LEC system with the objective of optimizing simultaneously the plant 
dynamic performance and minimizing the accumulated damage and/or damage rate in critical plant components. Figure 2. 1 shows 
a conceptual view of the LEC system. The high-level philosophy taken in this report is to develop an LEC system control approach 
that can easily be addended to a conventional performance controller design. 


FEED-FORWARD 



Figure 2.1. - Schematic diagram of Life Extending Control system and off-line optimizer. 
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The approach taken here is to design an aggressive performance controller to achieve a high level of dynamic performance. In 
the typical situation, with a linear or (linearizable) plant, the performance controller is designed using any of a number of linear 
controller design techniques (such as H oo , p synthesis, etc.). This assures stability in the inner loop with good dynamic performance 
as a starting point for the design of the highly nonlinear outer loop. This system of plant plus performance controller then becomes 
an augmented plant around which the life extending control loop is added. The essential elements of this outer loop are ( 1 ) a structural 
estimator that uses a set of plant outputs to estimate the load conditions (stress, temperature, or strain at the critical locations), (2) 
a damage model that uses these conditions to determine the rate of damage in the critical location, and (3) the damage controller. 
The damage model is a continuum time (as opposed to cycle) based representation of the damage so that it can be incorporated in 
the real time application. 

The objective is to reduce the damage rate and the accumulated damage at the critical points of the structure during transients 
where the time-dependent loads on the critical components can be controlled. This control action is usually indirectly applied by 
manipulating performance control inputs. The damage could derive from a variety of mechanisms such as microcracking wear, 
creep, fatigue, spalling, corrosion, and other mechanisms at one or more critical points. The time derivative of damage D indicates 
how the instantaneous load is affecting the structural components. The plant and remaining system dynamics in figure 2.1 are 
modeled by nonlinear differential equations which satisfy the local Lipschitz condition (Vidyasagar, 1 992) within the domain of the 
plant operating range. The structural model consists of the solution of structural dynamic equations representing the critical 
components under load conditions. This model may be a detailed representation of the structural dynamic behavior of critical plant 
elements or may be as simple as the isolated loads at the critical points determined from minor computations. A general structure 
of the plant and damage dynamics and their constraints is represented as follows: 

Plant Dynamics: 

^ 7 - = /(*(0.«(0); *( r o) = *o (2-1) 

at 

Performance Controller: 

x p (k + 1) = Ax p (k) + B[y CMD (k)~ y(*)] (2.2a) 

u p (k) = Cxp(k) + D[y CMD (k) - **)] (2.2b) 

Damage Controller: 

x d (k + \) = Ex d (k) + FIXk) (2.3a) 

u d (k) = Gx d {k)+HiXk) (2.3b) 

Structural Estimator 

q = fl(y) (2-4) 

Damage Model: 

D = f 2 (q) (2.5) 


3. The Reusable Rocket Engine 

This section contains details on the plant used in this study. The plant under control is a reusable bipropellant rocket engine 
shown schematically in figure 3.1. The propellants, namely, liquid hydrogen fuel and liquid oxygen, are individually pressurized 
by separate closed-cycle turbopumps. Pressurized cryogenic fuel and oxygen are pumped into two high-pressure prebumers which 
feed the respective turbines with fuel-rich hot gas. The fuel and oxidizer turbopump speeds and, hence, the propellant flow into the 
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main thrust chamber are controlled by the respective prebumer pressures. The exhaust from each turbine is injected into the main 
combustion chamber where it burns with the remaining oxidizer and is expanded through the rocket nozzle to generate thrust The 
oxygen flow into each of the two prebumers is independently controlled by the respective servo-controlled valves. The plant outputs 
of interest are the Oj/Hj mixture ratio the and main thrust chamber pressure which are closely related to the rocket engine 
performance in terms of specific impulse, thrust and combustion temperature. 



Figure 3.1.- Schematic diagram of reusable bipropeUant engine. 

A thermo-fluid-dynamic model of the rocket engine has been formulated for plant performance analysis and control systems 
synthesis (Ray and Dai, 1 995). Standard lumped parameter methods have been used to approximate the partial differential equations 
<V»«-rihtng mass, momentum, and energy conservation by a set of first-order differential equations. The plant model is constructed 
by '-""gal interconnection of the primary subsystem models such as main thrust chamber, prebumers, turbopumps, fuel and oxi r 

supply header, and fixed nozzle regeneration cooling. In this model, the plant has 18 state variables, two control inputs, and two 
outputs being controlled. The details of model development are presented in the appendix. 


4. Damage Modeling 

Damage modeling is a critically important aspect of life extending control. The damage model should have the following 
The model should be continuum based as opposed to cycle based for use in the control design process as well for 
implementation. Since the model is embedded in the life extending control loop it shouldbe as mathematically , and/or computanonaUy , 
simple as possible while representing the damage rate well enough to properly guide the actions of the controller. The implication 
of this is that the absolute level of the damage rate may not be so important as the form of the damage equation (or formulation). 
Further, computational simplicity becomes especially important when optimization is used in the design process. 

A wide variety of damage mechanisms are possible in the reusable rocket engine studied in this report. These mclude fatigue, 
spalling, high temperature creep, corrosion, and more. The objective of this report is to establish a viable design method for LEC 
systems containing extreme nonlinearities. Fatigue damage of the oxygen and hydrogen turbopump turbine blades is selected as the 
damage mechanism (and critical locations). As will be seen, this type of damage is extremely nonlinear and damage controller 
synthesis techniques which work for it will likely be adequate for other damage mechanisms. 

The fatigue damage model used in this study is a continuum-based analytical model (Lorenzo, 1994). This damage model offers 
two levels of treatment: ( 1 ) a local stress based model and (2) a more accurate strain-strain rate based model. Because of its simplicity 
and the easy availability of stress estimates, the stress based approach is chosen. For purposes of this study it is assumed that damage 
only occurs during tensile loading. Three specific damage rate models may be used to estimate damage rate D. 

For the case of zero mean tensile stress 
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-(1 +b)/b 


do . _ . 

— for o > 0,c increasing 


(4.1) 


where cris the instantaneous stress, and b and oj are material constants. When the mean stress is not zero, the damage rate is 
determined as 


v-l/i 


m=— 


Of — o„ 


(o-o m ) (1+b)/b — for a > a m > 0, o increasing 
dt 


( 4 . 2 ) 


where o m is the mean stress. Finally, under some conditions, it may be desirable to use a maximum damage rate equation, namely. 


2(1 + b)~ {l+b),b do , , 

ZKOm ax= — /or Of>O m > 0 . 


0 / _CJ 




( 4 . 3 ) 


This form provides an estimate of the instantaneous damage rate which is greater than the damage rate for any mean stress and is 
a conservative estimate for conditions when the LEC approach is to reduce peak stresses and mean stress is difficult to estimate. For 
the current application it will be seen that the damage mitigation is derived from reducing the mean stress on the turbine blades. 
Therefore, the mean stress damage rate equation (eq. (4.2)) is integrated to give the damage increment in one stress cycle as 


8 


eye 



^-1 tb 



( 4 . 4 ) 


where o a is the stress amplitude, o m is the mean stress, = 223.589 ksi is the fatigue strength coefficient, and b = -0.0858 is the 
fatigue strength exponent. It i s noted that the fatigue strength coefficient was adjusted by a factor of 0.82 to best match a more detailed 
damage model. The damage rate is calculated from the relation 


D = 


\-Ub 


av-a„ 


Q 

n 


(4.5) 


where £2 is the frequency of vibration of the blades in rad/sec. This model is used for both on-line damage estimation and in the 
optimization (off-line damage estimation). 


5. Design of the Dynamic Performance Controller 

The design of the inner loop performance controller is not the focus of this study. However, it must yield a well-behaved, stable 
closed-loop system compatible with the outer loop design process. It is designed to achieve aggressive dynamic response 
independent of damage considerations. The design procedure employed here uses the H ^ (or induced Z^ norm to Z^ norm) controller 
synthesis technique. This controller design method minimizes the worst case gain between the energy of the exogenous inputs and 
the energy of the regulated outputs of a generalized plant which is constructed below. Bamieh and Pearson ( 1 992) propose a solution 
to the induced Z^ norm controller synthesis problem for application to sampled-data systems. This design procedure has 
subsequently been incorporated as the function sdhfsyn in the MATLAB mutools toolbox (Balas et al., 1993). The performance 
controller needs to have very good low frequency disturbance rejection capabilities to prevent the damage controller output i/** 2 " 1 
from causing a long settling time in the plant outputs. 

Figure 5 . 1 shows the setup used for the synthesis of the induced Z^ norm controller for the rocket engine based on a plant model 
with two inputs (fuel prebumer oxidizer valve position and oxidizer prebumer oxidizer valve position) and two outputs (main thrust 
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chamber hot-gas pressure and Oj/Hj mixture ratio). The plant model is obtained by first linearizing the 18-state nonlinear model 
of the rocket engine at a combustion pressure of 2550 psi and an Oj/Hj ratio of 6.02. The bandwidth of the valves is assumed to be 
nonlimiting for this study. The pressure 2550 psi is chosen for linearization because the controller is required to operate in the range 
of 2100 to 3000 psi. After linearization, the 18-state linear model is reduced to a 13-state linear model for the controller design via 
Hankel model order reduction (Zhou, Doyle, and Glover, 1 996). A comparison of Bode plots reveals that reducing the 1 8-state model 
to 13 s t are* does not significantly alter the input-output characteristics of the original model. Since the induced Lj-norm controller 
synthesis procedure being used here requires a strictly proper generalized plant model, the problem of a nonzero D-matnx is 
circumvented by filtering the outputs of the controller by a first order filter with a very high frequency pole at 10 5 rad/sec, that is. 


W filter ( s ) - 


10 s 

s + 10 5 ’ 


(5.1) 



Figure 5.1. — Generalized plant. 


The frequency-dependent performance weight consists of two components: (1 ) W press , which penalizes the tracking error 

of combustion chamber pressure, and (2) which penalizes the tracking Oj/Hj error of the ratio. The frequency- 

dependent control signal weight W cont consists of two components: (1) W^which penalizes the fuel prebumer oxidizer valve 
position, and (2) W Q which penalizes the oxidizer prebumer oxidizer valve position. The objectives of the performance weights 
in this application are to keep steady-state error and overshoot/undershoot small while, at the same time, allowing a reasonably fast 
rise time. The objectives of these control signal weights are (1) prevention of large oscillations in the feedback control signal that 
may cause valve saturation and (2) reduction of valve wear and tear resulting from high-frequency movements. 

The parameters of both performance and control signal weights are initially selected based on the control system performance 
requirements and the knowledge of the plant dynamics; subsequently, the parameters are fine tuned based on the time-domain 
responses of the simulation experiments. For this design, the performance weights are 



(5.2) 



The control weight for both valves is 
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(5.4) 



Each of the two components of the frequency-dependent reference signal weight W^in figure 5.1 is chosen to be 


"V(5) = 


0.5 

5 + 0.5 


(5.5) 


A sampled-data controller, which is optimal in the induced Z^-norm sense, is designed using the generalized plant from figure 
5.1. As guaranteed by the design method employed, the controller has 21 states, which is the same as the number of states in the 
generalized plant model which consists of the reduced order plant model (13 states), the control signal filters (2 states), the 
performance weighting matrix (2 states), the reference signal weighting matrix (2 states), and the control signal weighting matrix 
(2 states). The controller provides acceptable reference signal tracking for the plant without using a large amount of control effort. 
It is found that reducing the order of the sampled-data controller from 21 states to 15 states causes no significant change in the 
controller dynamics from an input/output point of view. Therefore, this reduction causes no noticeable difference in the simulation 
results produced by the 21- and 15-state controllers. The 15-state controller is used in what follows. 


6. Damage Controller Design 

This section describes the design of the nonlinear damage-mitigating control loop. Here the output of a linear damage controller 
is added directly to the input of the plant, as shown in figure 2.1. The plant input is 

u(k) = Uff(k) + u ft> (k) + u dam (k) (6.1) 

where u^k) is a feedforward signal based on linear interpolation of steady-state values of the plant inputs, u^{k) is the output of the 
performance controller, and u ( j am (k) is the output of the damage controller. The damage controller is chosen to be a linear time- 
invariant discrete time structure. The procedure to be discussed below can be applied to damage controllers with a nonlinear 
structure; however, for simplicity, the procedure is demonstrated here using a linear damage controller. This section discusses a 
method which can be used to obtain the state space matrices (i.e., the A, B, C, and D matrices) of the Unear damage controller. 

The linear damage controller is designed by directly optimizing the elements of its A, B, C, and D matrices. To decrease the 
number of parameters to be optimized, the A matrix is constrained to be a diagonal matrix with distinct real elements. This is 
equivalent to constraining the damage controller to having unrepeated real eigenvalues. Repeated and/or complex poles can be 
included at the expense of computational complexity. For a damage controller with m inputs, p outputs, and n states, the number 
of parameters to be optimized is n (for the diagonal nxn A matrix) + nm (for the nxm B matrix) + pn (for the pxn C matrix) + pm 
(for the pxm D matrix) = n(l + m + p) + pm parameters. 

The parameters of the Unear dynamic filter are identified by minimizing a cost functional using nonUnear optimization. For each 
evaluation of the cost functional, a nominal computer simulation must be performed. The cost functional is evaluated by the 
simulation, and the simulation results are a function of the current damage controller chosen by the optimization routine. Since 
damage controllers designed using this method are directly based on the maneuver used in the optimization process, the maneuver 
should be chosen to be broadly representative of all plant operation. The resulting damage controller is then vaUdated by examining 
the results of various other typical maneuvers that the plant is expected to perform with this damage controller in the damage feedback 
loop. 

The simulation on which the design of the damage controller is based is a ramp-up of the main thrust chamber hot gas pressure 
from a level of 2700 to 3000 psi at a rate of 3000 psi/sec, followed by a steady state at the final 3000-psi pressure for 500 ms (see 
fig. 6.1). The 0 2 /H 2 mixture ratio for this simulation is to be kept at a constant value of 6.02. After each simulation is performed, 
data representing the results of the simulation are sent to the cost functional subroutine. These data consist of samples at every 
T = 0.002 sec of the chamber pressure, the 0 2 /H 2 mixture ratio, the damage rate in the 0 2 turbine blade, and the damage rate in the 
turbine blade. Since the duration of the simulation is 0.6 sec and each trajectory is sampled every T = 0.002 sec, there is a total 
of N = 300 samples sent to the cost functional subroutine for each of the four trajectories listed above. In addition, the value of 
accumulated damage for the 0 2 and H 2 turbines at time t = 0.6 sec is also used for the calculation of the value of the cost functional. 
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Figure 6.1. - Reference trajectory for chamber pressure. 


The cost functional includes the effects of both reference signal tracking (dynamic) performance and damage in the turbine 
blades; that is. 


jtot = J per f + jdam 

The performance part of the cost functional JP*** is the sum of penalties on 

(1) Tracking error of the main thrust chamber hot gas pressure (in units of psi) 


'press 


N - 1 

•lA 

k= 1 


press] 


~ y press {kT)-y T p j ess {kT) 

y%ss(kT) + 1.0 


( 6 . 2 ) 


(6.3) 


and 


J 5S = O 
J press 


ss 

press 


yprcssiM") y press ( NT) 


ypless(NT)- 


1.0 


(6.4) 


(2) Tracking error of the Oj/Hj mixture ratio 

N - 1 

J 0 2 /H2 =lQo 2 'H 2 g( kT ) ( 6 - 5 ^ 

*=1 


where 


-|2 


g(kT) = 


yo 2 /H 2 ( kT >- 602 

6.02 + 1.0 


0 


if yo 2 /H 2 ( kT ) s604 

otherwise 


( 6 . 6 ) 
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ho /h (NT) -6.02 f 
J ° 2 m 2 -Qo 2 /h 2 fi.02 + ,0 



(4) Accumulated damage in the H 2 turbine blades 

Jd„ 2 =Qd„ 2 [*>H 2 (NT)-Dh 2 (0T)\ (6.11) 

Both the pressure and 0^f^2 ra ^° components of the cost functional have extra weight on the error at the final sampling instant 
(i.e., the 14 th sample). Adjusting these extra weights is a means to control the steady-state behavior of the simulation. Increasing 
Qpress anc ^ or Qo^h 2 ten< * s to ^ ecrease settling time of the system. Also, since it is desirable to keep the O 2 /H 2 mixture ratio below 

a value of 6.04 during the transient, the Oj/Hj mixture ratio is penalized only if it exceeds 6.04 for samples 1 to N - 1. The final, 
N*, sample of the ra ^° * s penalized whether its value is above or below 6.04, since it is necessary for the O 2 /H 2 ratio to reach 

6.02 in the steady state. The factor of 1 .0 added in the denominator of equations (6.3), (6.4), (6.6), and (6.7) is a convenient way to 
combine the features of absolute and relative error and is often used in practice (Gill, Murray, and Wright, 1981). 

The accumulated damage and damage rate components of the cost functional do not contain an absolute value operator or 
squared terms because damage rate and accumulation are always positive. In the accumulated damage components (eqs. (6. 1 0) and 
(6. 1 1)), the initial accumulated damage is subtracted from the final damage at time NT = 0.6 sec to penalize the damage accumulated 
during the maneuver. The initial fatigue damage for both the 0 2 and the H 2 turbine blades is assumed to be D(0) = 0.1. 

Since the governing equations and the cost functional are nonlinear in nature, a nonlinear programming technique is used to 
identify the optimal parameters of the damage controller. Also, in order to evaluate the cost functional, a time consuming simulation 
must be performed. Therefore, a nonlinear programming technique known as Sequential Quadratic Programming (SQP) is 
employed; this technique has the reputation of being able to efficiently and successfully solve a wide range of nonlinear programming 
problems in which the evaluation of the cost functional is a computationally intensive procedure (Schittkowski, 1 985). A Sequential 
Quadratic Programming (SQP) Fortran Software package developed by Gill et al. (1986) at Stanford University called NPSOL is 
utilized to design the damage controller. 
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Interaction effects between the damage controller and the performance controller are minimized (1) by requiring a high level 
of dynamic performance through the cost functional for the nonlinear optimization of the damage controller, and (2) by the inherent 
frequency separation of the high frequency damage loop and the lower frequency performance loop. 

The following set of weights are found to produce an effective da m a g e controller: 

Q press = 21.0 Q boi =719.42 

Qpress - 10-6 Q Dh2 = 5.31 X 10 4 

Qo 2 l H 2 = 2.6x10 s Q Doi = 3.60x10 s 
Qo 2 ih 2 = 2-6 x 10 8 Q Dh2 = 2.66 x 10 7 

7. Simulation Results and Discussion 

The rocket engine is a two-input, two-output application (m = 2, p = 2), and the damage controller is designed using 15 states 
(n = 15). Therefore, the number of parameters to be optimized is 79. It is found that, after designing the 15-state damage controller, 
reducing the number of states to 5 via Hankel model order reduction does not significantly change the input/output characteristics 
of the controller. Therefore, the results that follow are created by using the 5-state reduced order damage controller. This result 
implies that it would be more efficient to directly optimize a controller with 5 states instead of 1 5 states. Unfortunately, it is not known 
a priori how to optimally choose the number of controller states. 

The damage controller is designed based on a transient which takes the chamber pressure from 2700 to 3000 psi (see figs. 7. 1 
to 7.6). Each plot displays the following cases: (1) no damage control (i.e., u(k) = Ujflt) + u p {k)) and (2) with damage control (i.e., 

«(*) = «/*) + « p (*) + «/*))• 

The chamber pressure trajectories for the two cases are compared in figure 7.1. The damage controller causes a slower rise time, 
a longer settling time, and less overshoot in the chamber pressure transient. The damage controller also causes the 0 2 /H 2 ratio to 
deviate farther from the desired value of 6.02 than the case with no damage control as seen in figure 7.2. However, the mixture ratio 
settles to 6.02 at steady state and remains within acceptable bounds throughout the duration of the simulation for both cases. 


c7 N 



Time (sec) 

Figure 7.1. - Main combustion chamber hot gas pressure (2700 to 3000 psi). 



Time (sec) 

Figure 7.2. - Oj/ H 2 mixture ratio (2700 to 3000 psi). 
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The damage rate and accumulation plots for the first 1 sec of the 2700- to 3000-psi simulation are shown in figures 7.3 to 7.6. 
Also, Table I summarizes the accumulated damage after this time interval for the two simulation cases (i.e., with and without damage 
control) for the two turbine blades. 



Time (sec) 

Figure 73. - Accumulated damage in Hj blade (2700 to 3000 psi). 



Time (sec) 

Figure 7.4. - Damage rate in H 2 blade (2700 to 3000 psi). 



Time (sec) 


Figure 73. - Accumulated damage in 0 2 blade (2700 to 3000 psi). 
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Figure 7.6. - Damage rate in 0 2 blade (2700 to 3000 psi). 


TABLE I. —ACCUMULATED DAMAGE (AFTER 1 sec) FOR 2700- to 


3000-psi SIMULATION 



Without damage control 

With damage control 

Ratio 

H 2 blades 

1.13X10' 5 

6.15X10" 6 

1.8 

0 2 blades 

1.21X10" 3 

3.45X10" 5 

35.1 


The loss of dynamic response of chamber pressure (fig. 7.1) and the modestly increased excursion in mixture ratio is the cost 
incurred for the improved damage performance. It is also observed that the slope of the accumulated damage (damage rate) at 
t = 1.0 sec for the H 2 turbine blade (fig. 7.3) indicates that there may be a relatively large steady-state damage rate for that turbine. 
If this is found to be the case for longer times, then the steady-state damage accumulation would far outweigh the transient da m age. 

The quality of the damage controller designed above is now tested on a transient maneuver which takes the chamber pressure 
from 2100 to 3000 psi at a rate of 3000 psi/sec (see figs. 7.7 to 7.12). This maneuver involves a larger pressure increase than the 
nominal maneuver used to design the damage controller, and, therefore, is expected to produce a larger amount of damage 
accumulation. 

A comparison of the chamber pressure trajectories with and without the damage controlleris shown in figure 7.7. As in the 2700- 
to 3000-psi case, the damage controller acts to “slow down” the transient as it approaches the final pressure of 3000 psi. Although 
the damage controller causes the Oj/Hj ratio to deviate from the desired value of 6.02 more than it did during the 2700- to 
3000-psi simulation, as seen in figure 7.8, it settles to 6.02 at steady state and remains within acceptable bounds throughout the 
simulation. The mixture ratio is important in this application as an indicator of chamber temperature (as well as propellant utilization) 
since the damage model does not contain temperature effects. Future implementations of the damage model can incorporate such 
effects. 



Figure 7.7 - Main combustion chamber hot gas pressure (2100 to 3000 psi). 
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Time (sec) 

Figure 7.8.- Oj/H 2 mixt ure ratio (2100 to 3000 psi). 

The damage rate and accumulation plots for the first 1 .2 sec of the 2100- to 3000-psi simulation are shown in figures 7.9 to 7. 1 2. 
Table II summarizes the accumulated damage for this transient 



Time (sec) 

Figure 7.9. - Accumulated damage in H 2 turbine blade (2100 to 3000 psi). 



Time (sec) 

Figure 7.10. - Damage rate in H 2 turbine blade (2100 to 3000 psi). 
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Time (sec) 


Figure 7.11. - Accumulated damage in 0 2 turbine blade (2100 to 3000 psi). 





Time (sec) 


Figure 7.12. Damage rate in 0 2 turbine blade (2100 to 3000 psi). 

TABLE H.— ACCUMULATED DAMAGE (AFTER 1.2 sec) FOR 2100- to 
3000-psi SIMULATION 


| | Without damage control 

With damage control 

Ratio 

H 2 blades 

2.46x1 O' 5 

9.61xl0" 6 

2.6 

O z blades 

2.48X10' 3 

7.01X10- 5 

35.4 


A deeper understanding of how the damage reduction is achieved may be obtained by observing the other state variables. The 
mechanism for damage reduction in this application is the reduction of mean stress on the turbine blades. This is achieved by the 
control by reduction of the peak value of turbine torque response (figs. 7.13 and 7.14). 



Time (sec) 


Figure 7.13. - Torque response in 0 2 turbine (2700 to 3000 psi). 
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Figure 7.14. - Mean stress in 0 2 turbine blade (2700 to 3000 psi). 


8. Summary and Conclusions 

The key concept of life extending control (LEC), as presented in this report, is to separate the design of the performance 
controller and the damage controller. A two-tier architecture has been proposed for the life extending control system which consists 
of a linear performance controller in the inner loop and a nonlinear damage controller in the outer loop. The high performance 
controller in the inner loop is designed using standard (linear) techniques (H^ or pi) to achieve an acceptable dynamic response for 
a reusable rocket engine which is similar to the Space Shuttle Main Engine (SSME). The combination of the rocket engine dynamics 
and the performance controller in the inner loop becomes the augmented plant for the design of the nonlinear damage controller (in 
the outer loop) which is the cascaded combination of a nonlinear characterization of fatigue damage rate in the turbine blades and 
a linear dynamic filter. The parameters of the filter are optimized to reduce the damage rate and accumulation at the critical points 
(i.e., fuel and oxidizer turbine blades) specifically under transient operations during which time the time-dependent load on the 
stressed structure is controllable. Benefits of this controller design approach are the following: (1) the performance controller can 
be designed by conventional (linear) techniques using commercially available software, (2) the effectiveness of the damage 
controller can be readily assessed relative to the reference design, and (3) when properly designed, the two- tier controller architecture 
can function over a broad range of transient requests (inputs) and not require an optimized feed-forward control sequence which is 
sensitive to plant modeling uncertainties and variations in the initial conditions. 

The damage controller designed in Section 6 reduces transient damage in the turbine blades of the reusable rocket engine by 
factors of 1.8 to 35 times as compared to when there is no damage controller in place. This reduction results in only a very small 
amount of degradation in the transient performance. 

Guaranteeing the stability of the closed-loop system is the single most important requirement of any control system design. For 
linear time-invariant systems stability can easily be determined by examining the eigenvalues of the A matrix of the closed-loop 
system. However, in general, proving the stability of nonlinear and/or time- varying systems is not very straightforward. In fact, for 
some complex systems it can be nearly impossible to analytically establish the stability of the system. Unfortunately, since fatigue 
damage processes contain severe nonlinearities, control systems containing a damage model in a feedback loop are nonlinear, and 
possibly time- varying as well. For the control system designed here, typical simulations were performed for which the closed-loop 
system was stable. Further, the apparent stability was increased as manifested by the chamber pressure response. However, for 
nonlinear systems, good performance and stability for a set of simulations does not guarantee that the system will be stable for a 
simulation not in that set For this reason a rigorous proof of the stability of damage-mitigating control systems should be pursued. 
Unfortunately, at this time, no proof of the stability of damage-mitigating control systems is available. A formal proof is a subject 
of future research. 
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APPENDIX 

THERMO-FLUID DYNAMIC MODELING OF THE REUSABLE 

ROCKET ENGINE* 

This chapter (Ray and Dai, 1 995) presents a nonlinear dynamic model of the thermal-fluid dynamics in a reusable rocket engine. 
The purpose of this model is to represent the overall dynamic performance and component interactions with sufficient accuracy for 
control synthesis and damage prediction. The governing equations used in the model are based on the fundamental principles of 
physics as well as on the experimental data under a variety of plant operating conditions. The model is formulated in the state-variable 
setting via nonlinear differential equations with time-invariant coefficients. 

The operating principles of the rocket engine under consideration are briefly described in Section A.l. Section A.2 presents 
the development of the nonlinear dynamic model equations using lumped parameter approximation. 

A.l Description of the Reusable Rocket Engine 

The reusable bipropellant rocket engine, under consideration in this report, is similar to the Space Shuttle Main Engine (SSME). 
Figure 3 . 1 in the mainbody of the report shows a functional diagram for operations and control of the rocket engine. The propellants, 
namely, liquid hydrogen and liquid oxygen, are individually pressurized by separate turbopumps. Pressurized liquid hydrogen and 
oxygen are pumped into individual high-pressure prebumers which feed the respective turbines with fuel-rich hot gas. The exhaust 
gas from each turbine is mixed in the common manifold and then injected into the main combustion chamber where it bums with 
the oxidizer to make most efficient use of the energy liberated by combustion. The oxygen flow into each of the two prebumers is 
independently controlled by the respective servo- valve while the valve position for oxygen flow into the main thrust chamber is held 
in a fixed position to derive maximum possible power from the engine. The plant outputs of interest are 0 2 /H 2 mixture ratio and 
combustor pressure which are closely related to the rocket engine performance in terms of thrust-to-weight ratio and engine 
efficiency. The liquid hydrogen is used as a regenerative coolant for the walls of the combustion chamber and thrust nozzle where 
structural integrity is endangered by the high temperature environment. The pressurized liquid fuel is circulated through the coolant 
jackets to absorb the heat transferred from the hot reaction gases to the thrust chamber and nozzle walls. 

A.2 Development of Plant Model Equations 

Standard lumped parameter approaches have been used to model the thermo-fluid dynamics of the engine in order to 
approximate the partial differential equations describing mass, momentum, and energy conservation by a set of first-order 
differential equations with time as the independent variable. The plant model is constructed via causal interconnection of the primary 
subsystem models such as the main thrust chamber, prebumers, turbopumps, valves, fuel and oxidizer supply headers, and 
regenerative cooling systems. The governing equations for the lumped parameter model of the plant dynamics are described in the 
following sections. In addition to the basic assumption of the lumped parameter approach, other pertinent assumptions are stated 
while describing the models of the individual subsystems. 

A A.l Fuel and Oxidizer Turbopump Subsystems 

The rocket engine has two sets of turbopumps, namely, low pressure and high pressure, for each of the two propellants. A 
simplified representation of the dynamic characteristics of the rocket engine is developed by lumping the low pressure and high 
pressure turbopumps into a single subsystem for each of the fuel and oxidizer propellants. On the oxidizer side, however, two pumps 
are modeled to obtain two sources of oxygen at different pressures. Model equations for the fuel and oxidizer turbopumps are given 
in Table A.l and Table A.2, respectively. 

Models of the hydraulic pump subsystems are derived based on the following assumptions: 

(a) The pump head which is proportional to the difference between static pressures at the suction and discharge is derived based 
on the assumptions of: (i) one-dimensional steady incompressible flow with negligible heat transfer; (ii) identical fluid velocities 
at the suction and discharge section of the pump; and (iii) no change in potential energy 
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Table A.l. — The fuel turbopump model equations 


Fuel Pump Model Equations 


Fuel Turbine Model Equations 


S pmp ” J 0 + S PMP (0) 

S PM P = ( ^TRB “ X pup ) f CpMPMI 
^ PMP = (1 CpMPw)(W HPBH + W opBfJ ) 
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C W I 
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Table A.2. — The oxidizer turbopump model equations 


Oxidizer Pump2 Model Equations 


Oxidizer Pump3 Model Equations 

w WO = < 1+ c OTRIw) 


Wop 3 = ( w hpbo + w opbo ) 

x ( w CMBO + w HPBO + w OPBO ) 


g OP3 =c OP3G w OP3 7 s OPMP 

g OP2 =c OP 2G w WO /s OPMP 


G OP3P =<I> OPMPP( G OP3) 

g OP2P = <i> OPMPP( G OP2) 


2 

G OP3D = C OP3P S OPMpG OP3P 

2 

g OP2D = c OP2P s OPMP g OP2P 
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p OP2PE = p OPMPS + G OP2D 
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g 0P2X = <p OP2X ( g OP2 ) 


2 
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2 
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Oxygen Turbine Model Equation 
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(b) The static performance of the pump is based on empirical characteristics (Rockwell, 1989) where the pump head AP pMp , 
power V pMP , and efficiency T| pMp are modeled as functions of the ratio of mass flow rate, Wp^ to pump speed S: 

APp M p oe S 2 <t»i(0); v pmp 06 S 2 O 2 (0); and T| pmp <*= S<J> 3 (0) (A.l) 

where © = W pMP /S, and the functions d>j , <J> 2 , and d> 3 are obtained from Rockwell (1989). Therefore, the outputs of the pump model, 
namely, pump discharge pressure, temperature, enthalpy, and torque, can be obtained from the pump characteristics and 
thermodynamic state relations. 

The governing equations for the turbine model are formulated under the following assumptions: 

(c) The working fluid in the turbine is a perfect gas and the expansion process in the turbine is adiabatic. For the ideal frictionless 
process, the following relationship holds: 


Tfc/TouUdeal = (Pin /Pout ( A - 2 ) 

where T is the absolute temperature, P is the pressure, the subscripts “in” and “out” respectively indicate the inlet and the outlet of 
the turbine, the subscript “ideal” stands for the idealized isentropic condition, and k is the ratio of the specific heats at constant 
pressure and temperature, which is assumed to be a constant within the operating range of turbine. 

(d) No loss of pressure and enthalpy occurs between the prebumer outlet and turbine inlet. That is, 

PpBR = ^TRB,in ’> ^PBR = ( A.3) 

(e) Flow through the turbine is assumed to be choked, and the kinetic energy of the fluid in the prebumer chamber is negligible 
such that the stagnation pressure and temperature, P* and T*, are respectively identical to the static prebumer pressure and 
temperature, P and T. Therefore, the mass flow rate through the turbine can be expressed as: 

C4S- = C-^L (A.4) 

V T TRB,in V T PBR * PBR 

where the coefficient C is calculated from the steady-state data. 

(f) The turbine efficiency and the output torque are obtained from the empirical characteristics of the turbine (Rockwell, 1989) 
as: 


TItRB “ 'nTRB^I 


ideal J 


Xtrb = ^TRB ideal ® 


' s ' 

V^HjdeaT ; 


(A.5a) 


(A.5b) 


where ideal (i.e., isentropic) enthalpy drop AH ideal is given as: 


^^ideal — C pi'll! 


T ' 

1 out, ideal _ ^ 

i rr> ~ ^p A m 

A in 




( p . , . \ 

r out, ideal 


k-1 

k 


(A.6) 


The outputs of the fuel and oxidizer turbine models, namely, turbine pressure, temperature, enthalpy, flow rate, and output 
torque are obtained from thermodynamic relations as delineated in Tables A.l and A.2, respectively. 
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The state variables in the fuel and oxidizer turbopump subsystems are respectively the shaft speeds S PMP and S OPMP . The power 
delivered by each turbine is equal to the sum of the power required by the propellant pump, and power losses in the bearings, gears, 
seals, and wear rings. Therefore, the dynamics of shaft speed in each turbopump are given in terms of the difference in torque as: 



where I is the moment of inertia and X indicates the torque. 


(A. 7) 


A.A.2 Prebumer Fuel and Oxidizer Supply Header Subsystems 

The model equations of the prebumer fuel and oxidizer supply header subsystems are listed in Table A.3. The equations of fuel 
flow to each prebumer are approximated to simplify the complexity of flow boundaries. The fuel flow to the two prebumers accounts 
for the mixture of the coolant flow from the primary nozzle cooling region and the primary nozzle bypass. The governing equations 
of the fuel flow through the prebumer header are derived under the following assumptions: 

(a) The prebumer fuel supply pressure Pppg is proportional to the fuel flow pressure at the main fuel valve. 

(b) Two coolant flows, namely, main chamber coolant flow (W^p) and primary nozzle coolant flow (W NOZF ), varies in 
proportion to the total fuel flow fWpyp). Since the coolant control valve position is held fixed, it is treated as fully open. Accordingly, 
the fixed nozzle bypass flow is obtained by subtracting the main chamber coolant flow and the nozzle coolant flow as: 


W CMBF = CcMBF W PMP 

(A. 8a) 

W NOZF = C NOZF W PMP 

(A.8b) 

W FNBP = W PMP “ W CMBF “ W NOZF 

(A.8c) 


By neglecting the dynamics due to fluid inertance in the flow passages, the above simplifications (a) and (b) reduce four differential 
pqn^rinnc of momentum conservation into four algebraic equations. This approximation only affects the model accuracy at high 
frequencies because of relatively small fluid inertance. 

(c) For one dimensional, incompressible uniform flow through a pipeline or valve and neglecting the body force, the friction 
pressure drop through a pipeline or valve is expressed as: 


D2 A 2 P 


C = f 


L. 1 
D 2A 2 


for pipeline 


(A.9a) 


D 2 A 2 Ra 


C' = f — — Ly, R A = = 
D 2pA 2 A A 


for valve 


(A.9b) 


The state variables of the preburner fuel and oxidizer supply headers are: 

• W HPBH and WjjpgQ (fuel mass flow rates into the fuel and oxidizer prebumers); 

• W opBH and W OPBO (oxidizer mass flow rates into the fuel and oxidizer prebumers). 

The derivatives of the above four state variables are obtained from conservation of linear momentum over a control volume of a 
pipeline, 


d_ 

dt 


(W) = 


P«, 



where p is the average fluid density and C f is the inverse of equivalent fluid inertance. 


(A. 10) 
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A.A3 Main Chamber Fuel Injector Subsystem 

The fuel injector mixes the two branches of fuel-rich exhaust hot-gas from the two turbines and a small amount of fuel from the 
combustion chamber coolant path. Model equations for the prebumers, main thrust chamber, and fuel injector are listed in Table A.4. 
The governing equations of the fuel injector subsystem are derived under the following assumptions: 

(a) The flow of an incompressible working fluid at a low Mach number (e.g., M< 0.3) is governed by the following relation 
(Blackburn et al„ 1960) by assuming that the subsonic velocities exist throughout the orifices: 


Q = uA = C d A^2(P in -P out )/p (volumetric flow rate) (A. 12a) 

W = Qp = C d ^2 (Pi,, -Pou t )/p (mass flow rate) (A. 1 2b) 

where p is the average density which is approximated as the gas density at the combustor. 

(b) The flow into the fuel injector manifold is the sum of two turbine exhaust flows, and W QTR , and main combustion 
chamber coolant flow W^p The manifold pressure Pp^.j is derived from Eq. (A. 12b) as: 

„ ( W TRB + W OTR + W CMBf) 2 , „ , A 1Tl 

PfINJ = ~2 + *CMB 

C dPCMB 


(c) The mixed gas temperature at the fuel injector manifold is obtained as a weighted average of the two turbine inlet 
temperatures, T PBR and T OPB , and the main chamber coolant flow temperature, Tq^. That is, T^j = C Q T pBR + CjTq PB + 
CjTq^bf w *>ere the coefficients, C d , C 0 , Cj, and C 2 are obtained from the steady-state data under normal operating conditions. 

A.A.4 Oxygen Control Valve Subsystem 

The nonlinearities of control valves are compensated by inducing the inverse characteristics of valves (Rockwell, 1 989) in the 
control signal such the valve command becomes proportional to the valve area under steady-state operations. The oxygen control 
valve subsystem model has two state variables, namely, fuel and oxidizer preburner valve rotary positions. The dynamics of each 
valve are represented by a first order lag as: 


-(Arfpv)= A rEPV ~ U - " V - < A - 15a > 

dt xppv 

-(Aropv)= Aropv ~ Uopv (A.i5b) 

dt t opv 

where Upp V and U opv are the commands to the oxygen control valves, and A^y and A ROpv are the effective areas of the oxidizer 
control valves, and T is the time constant of the respective valve. 

In solving the nonlinear optimal open loop control problem, the two commands Upp V and U OPV correspond to the decision 
variables in the nonlinear programming which are bounded above and below via specified constraints. 

A.A.5 Preburner and Combustion Subsystems 

The dynamic equations for the combustion process are developed by employing the principles of conservation of mass and 
energy with the following assumptions. 

(a) Conservation of momentum is satisfied by assuming that gas pressure and temperature in the combustion chamber are 
spatially unif orm although they are time-dependent, and the kinetic energy due to gas velocity in the chamber is negligible. This 
assumption is valid for a low-frequency dynamic representation, and precludes the process of high-frequency acoustic propagation. 
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Table A. 4. — Valve, preburner, combustion, and fixed nozzle model equations 


Oxygen Control Valve Model Equations 

Main Combustion Model Equations 

/ — \ 
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(b) One-dimensional unsteady flow in the combustion chamber is represented by a first order differential equation of the rate 
change of mixture gas density which is related to the mass flow into and out of the chamber via conservation of mass. 



w in -w out 

V CMB 


(A. 16) 


where is the volume of the combustion chamber. 

(c) The conservation of energy equation yields: 

A (C v VpT) = 2 + FW 0z - Q heat (A.17) 

where F is the energy liberated by per unit mass of oxygen from a macroscopic point of view of the chemical process where the 
reaction dynamics is assumed to be instantaneous. Qheat is the heat transfer rate from the control volume to the coolant channel wall. 

(d) Based on the thermodynamic relationship of the perfect gas law, the average gas temperature in the combustion chamber 
is given as: T CMB = Pcmb / (Pcmb r ) where R is A* characteristic gas constant Therefore, the derivative of the main chamber 
pressure is obtained by rewriting the energy Eq. (A.17) as: 

“ ( P CMB ) = ( W FINJ H FINJ + W CMBO H OP2E “ W NOZ H CMB + W CMBO F “ QcMBw)/( C V V CMb/ R ) (A.18) 

at 

(e) The flow through the nozzle throat is choked. 

The model equations of the prebumer and combustor are given in Table A.4. The six state variables in two prebumers and main 
combustion chamber are: 

• P pBR and R pBR : (Fuel prebumer chamber gas pressure and density); 

* Pqpb Rqpb : (Oxidizer prebumer chamber gas pressure and density); 

♦ Pcmb ^ dR OrfB : (Main thrust chamber hot gas pressure and density). 

The governing equations in prebumers are similar to those in the main chamber because of the similarity of the physical processes. 

A.A.6 Main Thrust Chamber/Fixed Nozzle Cooling Subsystems 

The basic relations governing the thrust chamber performance, such as specific impulse, combustion temperature and pressure, 
are calculated based on the thermodynamic principles of ideal rocket propulsion systems (Sutton, 1992). The following assumptions 
are used to derive the governing equations of heat transfer in the coolant channel wall. 

(a) The hot-gas velocity, pressure, temperature, and density are uniform across any cross-section normal to the nozzle axis. 

(b) No shock waves or discontinuities exist in the flow through the convergent-divergent nozzle, and the boundary layer effects 
are neglected. The energy equation applied across the nozzle throat and nozzle exit yields the exit temperature T c as a function of 
the throat temperature, T t , and exit Mach number M. 


T c =l 


V 


1 + — M 2 
2 




(A.19) 


where the exit Mach number M can be obtained as a function of the throat/exit pressure ratio, P/P e , and throat/exit area ratio, 
A,/A e , by combining the energy and continuity equations: 



(A.20) 
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In the simplified model of the main thrust chamber coolant channel subsystem, heat transfer rates and wall temperatures are 
derived using a lumped parameter model with two nodes. The model equations of the main chamber and nozzle regeneration cooling 
heat transfer subsystems are listed in Table A.5. The heat transfer process is characterized by three different mechanisms, namely, 
convective heat flux from the hot gas to hot- side of the coolant wall, the conductive heat flux through the wall from the hot-side to 
the cold-side, and the convective heat flux from the cold-side of the wall to the liquid coolant. 

(c) The conduction heat transfer rate is expressed in terms of a constant thermal conductivity of the coolant wall material and 
the temperature difference between the hot and cold sides as: 

q U c=(^)t w 2 -T w i ) (A.21) 

where k is the coefficient of thermal conductivity, and A is the area of heat transfer. 

(d) Convective heat transfer is associated with the mass transfer in a fluid boundary layer over a fixed wall. In. A.2, the rates 


of convective heat transfer and Q wf are given as: 



Qgw — h c A(T g - T w 2 j 

from the hot gas to hot-side wall 

(A.22a) 

Qwf ~ h c A(T wl — Tf ) 

from the cold-side wall to coolant 

(A.22b) 


Table A.5. — Main chamber and fixed nozzle regenerative cooling model equations 


Main Chamber Regenerative 
Cooling Model Equations 

Fixed Nozzle Regenerative 
Cooling Model Equations 

T cum = j 0 r* cum ( r ) dt + T cmw\ (° ) 

Tnozw = J 0 ^ NOZW ^^ + T nozw (0) 

Tcmwi = J 0 T cuwl ( *)& + T CUW2 (0) 

TnOZF = J 0 ^ NOZ1 S + ^vozf(^) 

Tcmbf = ^ + TcMBF (0) 



■'O CMBF 

QnOZW ^NOZWH^NOZ 7 NOZW ) 

QcMBW “ ^CMBWH^CMB ~~ ^ CMWl ) 

X 


xM 8 

CffOZFH “ CnOZFH + C NO zq\ T nozf ) 

QcMBWW ~ C C MBK (TcMWI ~~ Tcmw\) 

QnOZWF = ^NOZFH^NOZW ~ ^ NOZF ) 
1 I0J8 

QcMBWF = C € mbFh(^+ CcmbQ^CMBF ) 

X 

”nozf I 

X ( TcMWl ~~ TcmBF f 



* NOZW “ (QnOZW QnOZWF ) ' t'NOZWC 

^CMW2 ~ ( QcMBW ~ QcMBWW ) / CcMBWC 

TjvoZWF = ( Qnozwf ~^^nozf^p,h 2 

T CMWl = ( QcMBWW “ QcMBWF ) ! CcMBWC 
^CMBF = ( QcMBWF ^CMBF^PJi 2 

X ( Tpmpe “ T CMBF )) / CcMBFC 

X PM PE " ^ NOZF ))^ ^NOZFC 
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where h c is the convective heat transfer coefficient, T w2 and T wl represent hot-side and cold-side wall temperatures, respectively, 
and T f represents the bulk temperature of the liquid coolant. The convective heat transfer coefficient is described as a function of 
the fluid mass flow rate W and other system parameters at specified operating conditions using the following empirical correlation 
(Rockwell, 1989): 


h c oc w 0 ' 8 from the hot gas to hot-side wall (A.23a) 

h c oc (l + CTf )W 0 - 8 from the cold-side wall to coolant ( A.23b) 

For a thermal system composed of a material of density p, specific heat c p , and a constant volume V , the energy balance equation 
takes the following form: 


dT SWork 

P c p V ” = Qin (t) — Q ou t (0 + Qgen (1) ^ 


(A.25) 


where Q in or is the heat flux entering or exiting the control volume, Q gen is the rate of heat generated within the control volume, 

and 8Work/dt is the time derivative of the work done upon the control volume. 

Two wall temperatures at the two nodes on the hot and cold sides of the coolant channel wall, Tq v ^ w and Tq^^w, and 
hydrogen coolant temperature, T CMBW¥ , are the three state variables in the heat transfer model of the thrust chamber coolant channel. 
In reality, these state variables correspond to wall temperatures at the throat location where the heat flux is the highest and failure 
is most likely to occur. In contrast, the thrust chamber nozzle is relatively less prone to failure because of lower temperature. One 
lumped heat transfer node with two state variables is used to model the heat transfer through the nozzle coolant channel. The five 
state variables in the heat transfer model of the combustion and nozzle walls are: 

• and T CMW2 ^ the cold-side and hot-side temperatures of the combustor wall. 

• T cmbwf and T nozwf are coolant fluid temperatures in the combustor and nozzle. 

• T nozw is the average wall temperature of the nozzle. 

Derivatives of wall temperatures, Tq^j, Tq^^ and T N0ZW , are obtained via Eq. (A.25) as: 


— = ( t cmwi) - (Qcmbww ” Qcmbwf )/Ccmbwc 
dt 

— = (Tcmw2 ) = (Qcmbw “ Qcmbwf V^cmbwc 
at 


(A. 26a) 
(A.26b) 


■— = ( t cmwf) = [Qcmbwf + w cmbf c p,h 2 ( t pmpe ~ t cmbf)]/ c cmbfc (A.26c) 

The cold-side and hot-side temperatures, T CMW1 and T CMW2 , of the combustor wall are denoted as Tj and T 2 for brevity in the creep 
damage model in Chapter 3. The model steady-state results are presented in Table A.6. 
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Table A.6. — Steady-state model results 


Process Variables 

Symbol 

Unit 

100% 

Load 

(State Variables) 



Model 

Results 

Heat Balance 

Fuel turbopump shaft speed 

SPMP 

rad/sec 

3570.74 

3577.6 

Oxidizer turbopump shaft speed 

SOPMP 

rad/sec 

2917.49 

2849.4 

Main thrust chamber hot-gas pressure 

PCMB 

psi 

3000.0 

3006.0 

Main thrust chamber hot-gas density 

RCMB 

lb/in. 3 

1.3358d-04 

1.2673d-04 

Fuel prebumer hot-gas pressure 

PPBR 

psi 

4831.0 

4938.7 

Oxidizer prebumer hot-gas pressure 

POPB 

psi 

4854.09 

5003.5 

Fuel prebumer hot-gas density 

RPBR 

lb/in.^ 

4.7846d-04 

5.4478d-04 

Oxidizer prebumer hot-gas density 

ROPB 

lb/in.^ 

6.4924d-04 

6.7526d-04 

Fuel flow rate into the fuel prebumer 

WHPBH 

lb/sec 

82.1055 

78.18 

Fuel flow rate into the oxidizer prebumer 

WOPBH 

lb/sec 

76.1259 

67.78 

Oxidizer flow rate into the fuel prebumer 

WHPBO 

lb/sec 

38.5659 

35.1 

Oxidizer flow rate into the oxidizer prebumer 

WOPBO 

lb/sec 

20.665 

23.67 

Oxidizer flow rate into the thrust chamber 

WCMBO 

lb/sec 

809.656 

801.77 

Coolant side chamber wall temperature 

TCMW1 

°R 

1240.43 

/ 

Hot-gas side chamber wall temperature 

TCMW2 

°R 

1457.45 

/ 

Main thrust chamber coolant temperature 

TCMBF 

°R 

483.341 

469.1 

Coolant side nozzle wall temperature 

TNOZW 

°R 

1078.21 

1260.0 

Nozzle coolant temperature 

TNOZF 

°R 

433.145 

466.1 

Fuel prebumer oxygen flow valve position 

AFPV 

/ 

0.7813 

0.7812 

Oxidizer prebumer oxygen flow valve position 

AOPV 

/ 

0.6387 

0.6388 
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A3 Simulation of Transient Responses of the Rocket Engine 

In the thermo-fluid-dynamic model of the rocket engine derived above, the plant state vector consists of twenty state variables, 
two control inputs, and ten output variables as listed below: 


State Variables: 

Fuel Turbopump shaft speed 
Main thrust chamber hot-gas pressure 
Fuel prebumer oxygen flow valve position 
Fuel prebumer hot-gas pressure 
Fuel prebumer hot-gas density 
Fuel flow rate into the fuel prebumer 
Oxygen flow rate into the fuel prebumer 
Hot-side coolant wall temperature 
Cold-side coolant wall temperature 
Main thrust chamber coolant temperature 


Oxidizer Turbopump shaft speed 
Main thrust chamber hot-gas density 
Oxidizer prebumer oxygen flow valve position 
Oxidizer prebumer hot-gas pressure 
Oxidizer prebumer hot-gas density 
Fuel flow rate into the oxidizer prebumer 
Oxygen flow rate into the oxidizer prebumer 
Oxidizer flow rate into the main thrust chamber 
Nozzle cooling tube wall temperature 
Nozzle coolant temperature 


Control Inputs: 

Fuel prebumer oxidizer valve position Oxidizer prebumer oxidizer valve position 


Output Variables for Life Prediction and Plant Control: 

Main thrust chamber pressure mixture ratio 

Fuel turbopump shaft speed Oxidizer turbopump shaft speed 

Fuel turbopump torque Oxidizer turbopump torque 
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